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Abstract 

An unified approach for dealing with stability and aeroelastic response to time-dependent 
pressure pulses of swept wings in an incompressible flow is developed. To this end the indicial 
function concept in time and frequency domains, enabling one to derive the proper unsteady 
aerodynamic loads is used. Results regarding stability in the frequency and time domains, and 
subcritical aeroelastic response to arbitrary time-dependent external excitation obtained via the 
direct use of the unsteady aerodynamic derivatives for 3-D wings are supplied. Closed form 
expressions for unsteady aerodynamic derivatives using this unified approach have been derived 
and used to illustrate their application to flutter and aeroelastic response to blast and sonic-boom 
signatures. In this context, an original representation of the aeroelastic response in the phase - 
space was presented and pertinent conclusions on the implications of some basic parameters 
have been outlined. 
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Nomenclature 


Dimensionless elastic axis position measured from the midchord, positive aft 
c„ Chord length of wing, normal to the elastic axis, 

Lift-curve slope 

C{k\F{k),G{k) Theodorsen’s function and its real and imaginary counterparts, respectively 

/* , /„ Plunging and pitching deflection functions 

h , Plunging displacement and its amplitude, respectively 

H, , /t, Dimensionless unsteady aerodynamic coefficients 

L, , Mj Dimensionless unsteady aerodynamic complex coefficients 

i Imaginary unit, 4^ 

/,, , r„ Mass moment of inertia per unit length of wing and the dimensionless radius of 

gyration, (/ Imbl , respectively 
/ Wing semi-span measured along the mid-chord line 

Dimensionless aerodynamic lift, L^b„lmU^ 

Total lift and moment about the elastic axis per unit span of the swept wing 
L^,l^ Overpressure signature of the N-wave shock pulse and its dimensionless counterpart, 
L^b,^lmUl , respectively 

m,n Wing mass per unit length and wing/air mass ratio, mjKpbl , respectively 
Dimensionless aerodynamic moment, 

N Load Factor, \ + h"(g 

Pm, Pm Peak reflected pressure in excess and its dimensionless value , respectively 

r Shock pulse length factor 

s ,X Laplace transform variable and operator, respectively 

Static unbalance about the elastic axis and its dimensionless counterpart, 5^ Imb^ , 
respectively 

r ,To,T Time variables and dimensionless time, U„t!b ^ , respectively 

Freestream speed and its component normal to the elastic axis, respectively 
Dimensionless free-stream speed, 

X Coordinate parallel to freestream direction 

X Chordwise coordinate normal to the elastic axis 

y Coordinate perpendicular in the freestreem direction 

y Spanwise coordinate along the elastic axis 

w Downwash velocity 

z Transverse normal coordinate to the midplane of the wing 

Z Vertical displacement in z direction 

a Twist angle about the pitch axis and its amplitude, respectively 

Sr Tracer quantity 

Structural damping ratio in plunging, c^,/2/n&)^ and in pitching, , respectively 
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T] Dimensionless coordinate along the wing Span, >-// 

X , a Spanwise rate of change of twist and bending, respectively 
A Swept angle (positive for swept back) 

^ Dimensionless plunge coordinate, w!b„ 

p Air density 

Tp Dimensionless positive phase duration of the pulse, measured from the time of the arrival 
0(r),<I)(5)Wagner’s function in the time and Laplace domains, respectively 
(x),k„ Circular and reduced frequencies, cob^/U „ , respectively 

0)^,0)^ Uncoupled frequency in plunging, and pitching, , respectively 

o> Plunging-pitching frequency ratio, o)J(o„ 

Subscript 

( • )c Circulatory terms of lift and aerodynamic moment 
( • )„c Non-circulatory terms of lift and aerodynamic moment 
( • )„ Quantity normal to the elastic axis 
■ ) Quantity associated with the swept wing 

Superscript 

( . ) Variables in Laplace transformed space 

(•).(•) Derivatives with respect to the time / , and the dimensionless timer , respectively 


1. Introduction 

In this paper, the concept of the linear indicial''^ functions in the time and frequency domains 
is used to determine the associated unsteady aerodynamic derivatives for swept lifting surfaces. 
Such a treatment of the problem enables one to approach either the open/closed loop aeroelastic 
response in the subcritical flight speed regime to arbitrary time-dependent external excitations 
(such as e.g. gusts, airblasts due to explosions or sonic-booms'*’^), or the flutter instability of 
actively controlled/uncontrolled swept wings. In this paper both problems are addressed. 

As by-product of this analysis, a close form solution of the unsteady aerodynamic 
coefficients is obtained that is directly used in the stability and aeroelastic response problems. 

The unsteady aerodynamic lift and moment in incompressible flight speed regime are 
expressed for the swept aircraft wing in the time and frequency domains by using the Wagner 
and Theodorsen functions, respectively. For the response of dynamic systems it is only necessary 
to express the lift and moment via the indicial Wagner's function. For the approach of the flutter 
problem, the Theodorsen's function helps the conversion of the expressions of both the 
aerodynamic loads and the unsteady aerodynamic derivatives in the frequency domain®'^. 
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Herein the case of a 3-D lifting surface, including the plunging and pitching degrees of 
freedom is considered for flutter and response analyses. 

2. Preliminaries 

As shown in Ref. 8, for zero initial conditions, the unsteady aerodynamic loads can be 
converted from the time to the frequency domain via a Laplace transform. This results in the 
possibility of using the correspondence s ik„ as to convert the unsteady aerodynamic load from 
the time to the frequency domain, where s and k„ are the Laplace variable and the reduced 
frequency, respectively. The Laplace transform operator is defined as: 

^(- )=!“(• ( 1 ) 

In this sense, the Wagner's function <!>{t) is connected with the Theodorsen’s function c{k „ ) via a 

Laplace transform as: .^(O _ = <!>('<:„ ) , (2) 

ik„ ik„ Jo 

and vice-versa: <p{r) = X {c(A:„ )/ik „ } , Re[ik„ )>0 . (3) 

and having in view the correspondence s ik„ we can also write: 

) ,k„ "s r <t>{'^ydT = 0{s) . (4) 

Using this relationship, it is possible to obtain the full expression of unsteady aerodynamic 
coefficients in terms of the Theodorsen’s function C(k„) and its circulatory components 
F(/:„)and G(k„Y It is interesting to note that the reduced frequency parameter k„ for swept and 

for straight wings coincide: k„ = ^ and icot-ik^x. (5) 

fJ„ I/„cosA 

This implies that the indicial Wagner's function ^(t) remains invariant to any change of the 
sweep angle. 

3. Analytical Developments 

For swept wings, the total lift per unit span, can be expressed in the form: 

A (y - ^ )= A (y - ^ V- A 4cl A ^«<-2 ^ K A 0 • (6) 
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where the indices c and nc identify the various contributions associated with the circulatory and 
non-circulatory terms, respectively. Using similar notations, the total moment per unit span about 
the elastic axis is: 

aM ciya^AM „c2{y,ty^M (7) 

/^M„ca(yyf) being associated with the apparent moment of inertia^. Herein the lift is positive 
upward, while the moment is positive nose up. For the sake of convenience, herein the plunging 
coordinate is positive when is downward (see Fig. 1-3). Expressing the vertical displacement Z 
of a point on the center line of the wing as, (Fig.2): 

Z{x,y,t)=h{y,t)+xa{y,t), ( 8 ) 

where h = h[y,t), a = a(y,t) are the displacements in plunging and pitching, respectively, and 

assuming that the origin of the x axis coincides with the elastic center, the downwash velocity w 
normal to the lifting surface becomes: 

w(^x,y,t^= w(x,y,t^ = dZ/dt + U ^ . (9) 

The in-plane chordwise coordinate x normal to the elastic axis (see Fig. 2) can be expressed as: 

x=b„{ll2 -a„). ( 10 ) 

Consequently, using the dimensionless timer{=U jfb „ ) Eq. (9) becomes: 

w(x,y,T)=U — + a + —tanA + (l/2-a'{a' + b^^tanA\ ■ (11) 

Herein, b„ is the half-chord of the airfoil, U„ is the component of the flow speed, both normal to 
the elastic axis and ( )' = d{ )/dT . The quantity in Eq. (11) undescored by a solid line is usually 

discarded in the specialized literature'’'”, as being related to the wing camber effect. However, 
herein this effect will be taken into consideration. 

In the following sections, the unsteady aerodynamic loads in an incompressible flow can be 
obtained in time (3.1. a) and, with the use of the Laplace transform space in the frequency (3.1.b) 
domains. 
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3.1 Unsteady Aerodynamic Loads in Incompressible Flow 

3.1.a Time Domain 


The circulatory components of the lift expressed in terms of Wagner's indicial function ^(t) 
(referred also to as heredity function) obtained in time domain^ is: 

' h’ , d^h . f] ^ 


^ L, (y,T)= £ 0(t - To )| — + a' + tan A+ 


-a„ a +b, 


”^dTo 


tan A 


K, (12) 


As concern the aerodynamic non-circulatory components, using the dimensionless time these 
are expressed as: 

A ( F.t) = -| pUl [h’ - a„b„a] , ^ {y,x) = P^lb„a , 


\KcAy^'^)=~CLa,P^lbl tan A 


(5,+l)^ + 5,A + 5,— tanA 
hn 3y 


+ -a„C,„^pUlbli^nh 


(^r + UTT + ^rl^rtanA 


(13) 


Using the expression of the lift, Eq. (6), the equation for the moment, Eq. (7), can be cast as: 
a(yd) = -{y^ + a„yj„ ^L^(y,t)-a„b„^L„^.^{y,t) 

+ ( 1/2 AKciiy^t^AM „Ay<tyAM „Ay^t), ( 14 ) 

in which the last two non-circulatory components, using the dimensionless time, are expressed 
as: 


A (>’-''■)= P^lh] tan A + ic^„^ P^lbla^ tan A 


(5,+l)f + 5,A 


+ 5^ — tan A 
dy 


-\c,,pulb{ 


— + al ItanAl 
8 


(^r+l)7- + 5ri;rtanA 

b„ dy 


,M^{y,T)=~pC,„y„Ula'. 


(15) 

(16) 


Herein, the spanwise rates of change of bending and twist, cr and A, respectively, are expressed 
as a-dhldyaxid 'k = dajdy . In these equations, as well as in the following ones, the terms 
affected by the tracer 5^ are generated by the last term in the expression of the downwash 
velocity (Eq. (11)), (term underscored by a solid line). When these terms are discarded, being 
considered negligibly small, 5^ = 0 , otherwise 5^ - 1 . 

Replacement of Eqs. (12) and (13) in Eq. (6) and of Eqs. (13), (15) and (16) into Eq. (14), 
results in the unsteady lift and aerodynamic moment expressed in the time domain. Concerning 
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the circulatory part, Eq. (12), this can be explicitly determined by transforming this expression in 
the Laplace domain using the relationship between the Laplace transform of the Wagner and 

Theodorsen functions, namely C[-is)js = X (^(t)) = and afterwards inverting back the 

obtained expression in the temporal space. Alternatively, in order to ease the computations, the 
available approximate expressions for <^(t) and for can be used in the Laplace 

transformed space (see also Ref. 15 where the representation in terms of exponential polynomials 
and quasi polynomials was used to approximate the supersonic unsteady aerodynamic loads, 
representation that can be applied also for this case, as well). 

The expressions of lift and aerodynamic moment in the time domain, ^L^{y,T) and 
aM g{y,T) can be used to determine the subcritical aeroelastic response of swept wings. 
However, when the aeroelastic response of wings to time-dependent external excitations, such as 
gusts, sonic-boom, blast pulses, etc. is required, the unsteady aerodynamic loads in the time 
domain, and > have to be supplemented by the ones corresponding to above mentioned 
pulses. This will be considered in the next developments and an illustration of the capabilities 
provided by this unified method will be given. 

3.1.b Frequency Domain 

Upon replacing s^ik„ in the Eqs. (6) and (14) converted in Laplace transformed space; 

using the relationship between Laplace transform of Wagner and Theodorsen’s functions (Eq. 
(2)); representing the time dependence of displacement quantities as: 

a(y.T)=/„(y)a(T,^„) = /a(y)ao<;'*"\ h(y,r) = ff,[y)h{T,k„) = (17 a) 

and expressing: 

^L„(y,k„,T)=J^(y,ky^"\ (17 b) 

the equations for the unsteady lift and moment amplitudes can be expressed in the frequency 
domain (Appendix A). These expressions in the frequency domain can be used in the flutter 
analysis of swept aircraft wings. These coincide with the ones obtained differently in Ref. 16. 

In this analysis (y) and /y, (y) are chosen to be the decoupled eigenmodes in plunging and twist 

of the structure, and are determined as to fulfill identically the boundary conditions. These are 
expressed as: 
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fh (y)= iv = yli)=c\ +sin^i 

cosh "f" cos jljj 


( 18 ) 


/a(j)=/^a(^ = y//)=Qsin;32^, (19) 

where for the first bending and torsion we have p, = 0.5969n and P 2 = • The constants C, 

and Cj are chosen as to normalize fh{y) and /„(y), and so to get the unitary maximum 

deflection at the wing tip. The uncoupled first bending and torsion mode shapes were needed for 
the evaluations of the terms in the Eqs. (A. 1) - (A. 8) and are shown in Ref. 17. 


3.2 Unsteady Aerodynamic Derivatives in the Frequency Domain 

In this section close form solution of unsteady aerodynamic coefficients have been derived 
and their use in the process of unifying the aeroelastic formulations for flutter and response 
analysis of swept aircraft wings have been emphasized. 

A careful inspection of equations for lift and moment expressed in the time domain, Eqs. (7) 
- (8), suggests the following representations for the lift and aerodynamic moment: 


, ^ 2 ^ 7 7 /, 7/4 -^ + //5a' + //(, ^ 1 
2 b„ b„ 


. Ma (y. ^ P^n 2^„^ [ A 7“ + ^2«' + ^ + A 7" 


Herein, Hi, Ai denote the dimensionless unsteady aerodynamic coefficients, and where k„ has 
been included as to render the quantities in brackets nondimensional. In a simplified context, 
such a mixed form of the lift and moment was used in Refs. 18 and 19. The unsteady 
aerodynamic derivatives for swept wings are obtainable from the previous equations of lift and 
aerodynamic moment, by assuming harmonic time dependence of displacements quantities. In 
such a way, the frequency domain counterpart of Eqs. (28) - (29), expressed in compact form, 
becomes: 


A/ 1 cx,qM 2 ■> (22) 

Herein, the unsteady aerodynamic complex coefficients Lj and Mj can be expressed in terms of 
unsteady aerodynamic derivatives as: 

Lj=iH,+H^, L2=iH2+H2, M,=iA, + A^, M2=iA2 + A2- (23) 


L, k„ ) = pU^„ klb,^ L, + a^l2 , ^ (y, k„ ) = pf/ 
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where, for the sake of convenience, these are written as: 

H,=H,,H 2 =H 2 ,Hs= {Ms - Hs ), (24 a) 

and A, = A,,A,=A 2 ,A,= {A, - A, ). A, ={A,~ A,). (24 b) 

The close form solution for the unsteady aerodynamic derivatives in the frequency domain for 
swept wings will be obtained from Eqs. (22), expressed in terms of Wagner’s function 0(ik„) . 

Separating the real and the imaginary parts of these expressions, the unsteady aerodynamic 
derivatives result under the form displayed in Appendix A. These include also the spanwise rates 
of change of bending and twist, a and A, in terms of derivative in space of the shape functions, 
which are associated with the sweep effect. For straight wings these terms become immaterial. 
Notice that, the unsteady aerodynamic loads coincide with the ones obtained differently in Ref. 
16, and when specialized for 5^=0 , with the ones of Ref. 1. The direct use of the aerodynamic 
derivatives in the Eqs. (20) and (21) will help us to solve the eigenvalue problem associated with 
determination of the flutter instability. On the other hand, their Laplace transformed counterpart 
is helpful in the aeroelastic response analysis, wherein the equations of motion are converted to 
an algebraic systems of equations in the Laplace space. Notice that the flutter analysis can be 
conducted also in the Laplace space domain. In this case, classical methods such as U-g and p-k 
methods can be used in the Laplace transformed space. 


3.3 Aeroelastic Response of a Swept Wing to Blast and Sonic-Boom 
Pressure Pulses. Flutter Instability Derived from the Response. 


An application on the flutter instability and the aeroelastic response of a swept wing to blast 
and sonic-boom pressure pulses will be given in the next developments. The aeroelastic 
governing system of equations of a swept wing featuring plunging and twisting degrees of 
freedom to blast pressure signatures, expressed in dimensionless form, can be cast as: 


El 




dy 


w „ B 6 


,-4 




d^d 




-M^{y,t)=0. 


(25) 

(26) 


For the cantilevered wing, the related boundary conditions are: 
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dy 


— 0 and 


j=o 


( d^w{y,t) _ 9 

df df By 


= 0 (27) 


y=; 


Use of shapes functions, given by Eqs. (18) and (19), the aeroelastic governing equations in 
dimensionless form become: 

Xa^'i^)+ + /„( t )= /*( t ), (28) 

feAa>lT)+0'(O+(2C„/V„)9(T)+0/K,f-m„(T) = O, (29) 

The nondimensional parameters appearing in the preceding equations are displayed in Appendix 
B. Moreover, all these quantities are based on sections normal to the reference axis. 

The dimensionless sonic-boom overpressure signature of the N-wave shock pulse, can be 
expressed as follow; 


lbir)=H{r)^ 


m 


f \ 

1 -^ 




1 - 


(30) 


Herein, the Heaviside step function h{x) has been introduced in order to describe the typical 
pressure time-history for sonic-boom loads; denotes the dimensionless peak reflected 
pressure in excess of the ambient one (see Refs. 4,5,20 and the references therein); Zp denotes the 
positive phase duration of the pulse measured from the time of impact of the structure; r denotes 
the shock pulse length factor. For r = / the N-shaped pulse degenerates into a triangular pulse 
that corresponds to an explosive pulse (Fig. 4.a), and for r-2 a symmetric N-shaped pulse is 
obtained. A depiction of 4/p„, vs. time is displayed in Fig. 4.b. 

Equations (28) and (29) can be converted to the Laplace transformed space and solved for the 
unknowns, ^ {= X (<^)) and 6 {= X (0)); inverted back in time domain one obtain the plunging 

and pitching time-histories and the load factor time-history due to the sonic-boom pressure pulse, 
^(t) = '^(|(s)) and 6{r) = X ‘ (0(5)), respectively. Notice that, in this unified approach, when 


Eqs. (28) and (29) are converted in the Laplace transformed space, the expressions of 
aerodynamic derivatives expressed in the Laplace space can be used directly. 

Moreover, when the dynamic response of the actively controlled lifting surface is analyzed, 
also the feedback control forces and moments, that are time dependent, have to be included in 
Eqs. (25) and (26). 
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However, when the concept of the boundary feedback control^'’^'^ is used, the boundary 
conditions at the wing tip, Eqs. (27) have to be completed by the related control forces or 
moments. 


4. Results and Discussions 


Using the idea developed in Ref. 24, a modified strip theory has been developed as to address 
the problem of the aeroelastic response (for open/closed loop aeroelastic systems), by capturing 
also the 3-D effects. For swept wings, the local lift-curve slope for sections normal to the 

elastic axis are obtained from the aerodynamics of swept wings^'*, and are expressed as 

=Qa/cosA where =. 


AR cos A 


cos AC^„ jn + AR^\ + (Cj.„ cos A/(;c4/?))^ 
in which also the correction for the aspect ratio A/? has been included 

As a by-product of this analysis, a general close form solution for the unsteady aerodynamic 
derivative for swept wing are also included. Moreover, in the case in which also the 
compressibility effects are considered, the respective indicial functions for the compressible 
subsonic, supersonic and hypersonic flight speed regimes, have to be modified as to include also 
these effects. 

In general, the corrective terms identified by the tracer do not modify the trend of 


coefficients. For small aspect ratio wings, the maximum influence of the corrective term is 
present in the first plunging coefficients . Usually, for all coefficients, the effect of these 
terms becomes larger for higher sweep angles and for larger values of the parameter a „ . 
Whereas, for smaller aspect ratio wings, the effect of the camber becomes significant, and as a 
result it should be included, for higher aspect ratios, this effect becomes negligibly small. 

The aspect ratio AR parameter plays an important role on the coefficients H^, A^ and A ^ , and 
because these are correlated to the plunging displacement, the same role will be played on the 
aerodynamic loads. 

Notice that, for the circulatory components of the Theodorsen’s function assume the 


values F{k)-^ 1/2 and G{k)-^ 0 , and the corresponding unsteady aerodynamic derivatives can be 


determined, this being in agreement with the steady state solution for lift and aerodynamic 
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moment. In these developments, all the terms, including the aerodynamic ones associated with 
h and d , have been retained. Usually, these terms are neglected but, due to the presence of high 
frequency components in the blast pressure terms, their effect can be significant. 

As a result, the coefficients and also subsist. Whereas the aerodynamic 

coefficients /?, and Aj are the principal uncoupled aerodynamic damping coefficients in 
plunging and torsion, respectively, H 2 and A, are the coupled damping coefficients. As 
concerns, the depiction of H, and A, versus In/k , this representation enables one to get an idea 
of the variation of the respective quantity with that of the normal freestream speed . It should 
be mentioned that the expressions of the lift and aerodynamic moment in the frequency domain 
obtained using the correction of Ref. 16 coincide with the ones obtained here via indicial 
function approach. These coefficients are directly applied to the flutter analysis and the results 
are displayed in Fig. 5. The comparison concerns the flutter prediction via indicial function of a 
cantilever metallic swept wing of A = 4 and A = 30°. The results reveal that the predictions 
provided by the present approach, i.e. U p =235 ft js and oj/r = 386 raJ/s are in excellent 
agreement with those displayed in Ref. 1. In this graph the dotted lines represent the results 
based on the assumptions in Ref 1 whereas, the solid line represent those from the present 
analysis including also the camber effect. The critical value of the flutter speed is obtained here 
from the eigenvalue analysis of the homogeneous system of equation and at the same time from 
the response analysis. The way to determine the flutter speed from response becomes clear from 
the phase plane portrait (see Figs. 11-14), in which, for a certain restricted range of the flight 
speed (in the vicinity of the flutter speed), a periodic response with constant amplitude is 
experienced. 

On the same figure, determined from the response time-histories, there is also indicated the range 
in which the flutter instability occurs. The flutter predictions based on both methods show an 
excellent agreement. 

The graphs depicting the aeroelastic response time-history to blast pulses (i.e. explosive and 
sonic -boom blasts) are displayed in Figs. 6 - 14. In each of these, the corresponding type of 
blast/gust signature was indicated as an inset. In addition, the parameters in use for the 
simulations, unless otherwise specified, are chosen as; F„=l ; // = 10; m = 0.5; Tq = 0.5;x« = 0. 125; 
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Cti = Ca = 0;a„=-0.2) and (pm=l;Xp =15”;r=l or 2) for the structural and the load components, 
respectively. 

In Fig. 6 predictions of the aeroelastic response of a straight/swept wings to a gust load 
(indicated in the inset of the figure) are presented. Within this plot both the exact and selected 
approximate expressions of the Theodorsen’s function (Ref. 1) are used. As is readily seen, the 
differences occurring as a result of these approximations in the plunging time-history are 
indiscernible, which allows us to infer about the high accuracy of the approximations involving 
the expression of C(j). Nevertheless, for more complex wing configurations these approximate 
expressions can be highly useful towards determining the aeroelastic response. 

The graphs in Fig. 7 supply the dimensionless plunging {^ = wlb„) and pitching {o) 
displacements, and the load factor {N = \ + w"! g , where g is the acceleration of gravity) time- 
history aeroelastic responses to blast pressure pulse. It becomes apparent that an increase of the 
wing sweep angle, results in a decrease of the severity of the pulse signature. 

Moreover, the plunging-pitching coupling helps to reduce the amplitude of the aeroelastic 
response”. The load factor N has its maximum for t = 0, when the first impulse due the blast 
load occurs. 

Figures 8 highlights the effect of the speed parameter {=Ujb„(o„) on the swept aircraft 
wing (A = 15”) subjected to blast pulses. It becomes apparent that the amplitude of the response 
time-history increases with the increase of V„ . Moreover, in a certain range of speeds, as time 
unfolds, a decay of the amplitude is experienced, which reflects the fact that in this case the 
subcritical response is involved. However, for the dimensionless speed parameter V„ =2.1, the 
response becomes unbounded implying that the occurrence of the flutter instability is impending. 
The effect of the mass parameter n (= m/npbl) for swept wing of A = 15“ is indicated in Figs. 9. 
The increase of the mass ratio and consequently of the aspect ratio AR, (see Appendix A) results 
in the increase of the plunging and pitching displacement amplitudes. Moreover, for higher mass 
ratios and AR, the motion damps out at larger times. 

It should be noticed that the response to sonic-boom pressure pulses involves two different 
regimes (see Fig. 10); namely; one for which 0 < t < 30” that corresponds to the forced motion, 
and the other one to t > 30” belonging to the free motion. The jump in the time-history of N is 
due to the discontinuity in the load occurring at t= 30”. This jump doesn’t appear for explosive 
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pressure pulses, where r = 1 (see Figs. 7-9). Moreover, also in this case, when increasing the 
sweep angle, the effect of the blast signature becomes less severe. In Fig. 1 1 a three-dimensional 
plot depicting the dimensionless plunging ^ and pitching 9 deflection time-histories of a swept 
aircraft wing (A = 15°) to blast pressure signatures vs. the normalized spanwise coordinate y and 
the dimensionless time T has been supplied for two different values of the dimensionless speed 
parameter V„ , 1.00 and 2.5. 

The evolution of the aeroelastic system can be graphically illustrated by examining its 
motion in the phase - space, rather than in the real space, and recognizing that the trajectory 
depicted in this space represents the complete time history of the system. Corresponding to the 
flutter speed, that coincides with that otained from the eigenvalue analysis, the trajectory of 
motion describes an orbit with constant amplitude, the so called center. For V' < V/r as time 
unfolds, a decay of the amplitude is experienced, which reflects the fact that in this case a 
subcritical response is involved (stable focal point), while for V > Vf the response becomes 
unbounded implying that the occurrence of the flutter instability is impending (unstable focal 
point). 

Figure 12 highlights a three dimensional phase space portrait vs. and N ) of the 
plunging time-history response to blast load of a swept aircraft wing (A = 15°) for selected 
values of the speed parameter. As it is shown, the response becomes unbounded for = 2.1, 

representing the critical speed in which the periodic solution has been obtained and implying that 
the occurrence of the flutter instability is impending. To avoid its occurrence there are two 
possibilities, namely, including a passive/active control methodology or acting on the sweep 
angle A. In the latter case the idea is to use the capabilities of the aircraft featuring variable 
sweep angle (e.g. F-14 Tomcat) as to reduce the oscillations and at the same time expand the 
flight envelope. 

Figure 13 shows the phase plane portrait and the relative 3-D plots (vs. the load factor N) for 
selected values of the sweep angle A. With the increase of the sweep angle, the motion damps 
out at smaller times, the amplitude of the response is lower and the load factor is lower as well. 

In Figs. 14 a 3-D pictorial view of the motions in plunging and pitching around the flutter 
instability boundary vs. the variation of the sweep angle are displayed. From these plots a 
complete view of how the maximum maximorum of the admissible values of the amplitudes as 
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function of the sweep angle are evolving. It clearly appears that acting on an aircraft featuring 
swept wing angle variability, the subcritical aeroelastic response can be controlled in terms of the 
plunging and pitching displacement amplitudes. 

5. Conclusions 

An unified treatment of the aeroelasticity of swept lifting surfaces in time and frequency 
domains has been presented, and the usefulness, in this context, of the aerodynamic indicial 
functions concept was emphasized. The time domain representation is essential towards 
determination of the dynamic aeroelastic response to time dependent external loads, and in the 
case of the application of a feedback control methodology, of the dynamic aeroelastic response 
to both external time dependent loads and control inputs. The aeroelastic response have been 
represented in both, the classical way, by displaying the time-histories of plunging - pitching and 
load factor, and in an original phase - space context, that provides full information about the 
behavior of the aeroelastic system. The frequency domain representation is essential towards 
determination of the flutter instability. It was also shown how to capture the flutter instability 
also from the response time-histories. In this sense, the inclusive approach used here will be 
helpful toward the validation of the aeroelastic model via the comparison of the theoretical and 
flight vibration test results^^. 

Applications assessing the versatility of the methodology presented here toward the approach 
of both the subcritical aeroelastic responses and flutter instability for 3-D swept aircraft wing 
have been presented. The concept of the stability boundary and control of it via the use of the 
variable-sweep-wing-geometry have been illustrated. The unified formulation presented in this 
work can be extended as to approach it in various flight speed regimes^’’^*^ and to include also a 
feedback control mechanism. Moreover, the methodology developed in this paper can assist 
determination of the critical flutter speed from the investigation in flight or in wind tunnel tests 
of the aeroelastic response of aircraft wings to pulse loading. 
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Appendix A 


Unsteady Aerodynamic Derivatives: 
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where the aspect ratio AR is defined as AR-lJb„cos.^K and the spanwise dimensionless 
coordinate rj-y/ 1 have been introduced. For A = 0 the expressions of aerodynamic 
coefficients corresponding to straight wings are obtained. 



Appendix B 

Nondimensional parameters for flutter and response analyses: 
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Figure Captions 

Fig.l Nonuniform swept wing. 

Fig. 2 Airfoil section. 

Fig. 3 3-D view of a swept wing. 

Fig. 4 (a) Sonic -Boom and (b) Triangular Blast Pressure Pulses. 

Fig. 5 Flutter calculation via U-g method (Ref. 1). 

Fig. 6 Predictions of the aeroelastic responses of a straight/swept wings to a gust load, based on 
exact and selected expressions of the Theodorsen’s function. 

Fig. 7 Influence of angle of sweep A on the aeroelastic response to a blast pressure pulse at the 
wing tip. 

Fig. 8 Influence of the speed parameter on the response of a swept aircraft wing (A = 15°) to 
blast pulses at the wing tip. 

Fig. 9 Influence of the mass parameter |i on the response of a swept aircraft wing (A = 15°) to 
blast pulses at the wing tip. 

Fig. 10 Influence of the sweep angle A on the response to sonic-boom pulses at the wing tip. 

Fig. 11 Three-dimensional plots depicting the dimensionless plunging ^ and pitching 
0 deflection time-history of a swept aircraft wing (A = 15°) to blast pressure signature, vs. the 
normalized spanwise coordinate y and the dimensionless time T. 

Fig. 12 Three-dimensional phase-space portrait depicting the dimensionless plunging deflection 
time-history of a swept aircraft wing (A = 15°) to blast pressure signature, vs. the load factor N 
for selected values of the speed parameter. 

Fig. 13 Phase-plane phase-space portraits depicting the dimensionless 

plunging deflection time-history of a swept aircraft wing to blast pressure signature, for selected 
values of the sweep angle. 

Fig. 14 Cone of stability orbits depicting the envelope of the upper bound values of the 
dimensionless plunging and pitching deflection time-histories vs. A to blast pressure signatures. 
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